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It is thought that the current phase of accelerated expansion of the Universe is driven by a 
hypothetical fluid with negative pressure, dubbed dark energy. This fluid, of unknown nature, would 
correspond to approximately 70% of the energy content of the Universe. The cosmological constant, 
which may be associated with the zero point energy of the quantum fields, is the strongest candidate 
for dark energy, despite the problems that it brings to the cosmology-particle physics interface. The 
main alternative to the vacuum energy is a dynamical scalar held. From the observational viewpoint 
it is extremely hard, if not impossible, to decide which of these scenarios is correct, if one of them is. 

The data shows only that the Universe is expanding at an accelerated rate but does not reveal which 
causes this acceleration, dark energy or something else. However, although unknown, dark energy is 
not immune to the laws of physics. Such an exotic fluid must satisfy the bounds imposed by the laws 
of thermodynamics which have a strong experimental basis. Here we investigate the limits imposed 
by thermodynamics to a dark energy fluid. We obtain the heat capacities and the compressibilities 
for a dark energy fluid. These thermodynamical variables are easily accessible experimentally for 
any terrestrial fluid. The thermal and mechanical stabilities require these quantities to be positive. 

We show that such requirements forbid the existence of a cosmic fluid with negative constant EoS 
parameter which excludes vacuum energy as a candidate to explain the cosmic acceleration. We also 
show that the current observational data from SN la, BAO and H(z) are in conflict with the physical 
constraints that a general dark energy fluid with a time-dependent EoS parameter must obey which 
can be interpreted as an evidence against the dark energy hypothesis. Although our result excludes 
the vacuum energy, a geometrical cosmological term as originally introduced by Einstein in the held 
equations remains untouched. 

PACS numbers: 98.80.Es, 98.80.-k, 98.80.Jk 


Since the standard cold dark matter model (SCDM) 
was discarded by the Type la Supernovae observations, 
which point to a present day accelerated Universe [J, 2), 
physicists have been challenged to find a model that 
agrees with the data and, at the same time, rely on a 
solid theoretical basis. This task is not easy. Usually, 
physicists start with the simplest model. Thus, the first 
attempt was to reintroduce the cosmological constant A 
to Einstein’s field equations. A positive A term acts in 
the motion equations as a constant repulsive force and, 
therefore, can speed up the Universe at large scales. The 
so-called ACDM model is able to explain most of the 
current observational data and has a strong theoretical 
appeal since it may be linked with the zero point en¬ 
ergy of all quantum fields filling the Universe. Physicists 
would have given their verdict in favor of the cosmolog¬ 
ical constant if it did not suffer from a serious problem: 
the value of the vacuum energy density obtained by ob¬ 
servations differs from the value provided by quantum 
field theory by at least 60 orders of magnitude Q-Q- It 
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is very difficult to handle this huge discrepancy. If we 
think in terms of a net cosmological constant as the sum 
of a bare geometrical A term with the quantum vacuum 
energy density to explain such small value, this will gen¬ 
erate a fine-tuning problem: the absolute value of the 
geometrical and matter contributions to the net cosmo¬ 
logical constant must be extremely close. Also, symmetry 
arguments are not enough to explain the small value of 
the vacuum energy observed today. The lack of a rea¬ 
sonable explanation for the cosmological constant prob¬ 
lem has led physicists to explore other routes to explain 
the observations. Scalar fields (quintessence) stand out 
among the alternatives to the cosmological constant since 
they provide a link between particle physics and cosmol¬ 
ogy. The energy density of such a form of scalar matter 
must evolve with time but should mimic a cosmologi¬ 
cal constant to be compatible with the data. However, 
the majority of scalar fields that adjust the data have no 
foundations on particle physics and are somewhat arti¬ 
ficial. In fact, vacuum energy and quintessence are not 
the only possibilities to explain the cosmic acceleration. 
Other fluids with negative pressure called generally dark 
energy (DE) jf| and modifications of gravity theory Q 
are also in the game. 

By assuming that general relativity is the right the¬ 
ory of gravitation, the DE pressure must be sufficiently 
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negative to make the sum p -f 3p negative in order to 
produce an accelerated expansion of the Universe. Dark 
energy is frequently characterized by the equation of 
state (EoS) parameter (the ratio between its pressure 
and its energy density) w = p/pc 1 2 which can be con¬ 
stant or time-dependent. Such a phenomenological ap¬ 
proach enccompasses vacuum energy (w = —1), scalar 
fields (—1 < w < 1) 0, phantom fields (w < —1) @ 
and many other forms of exotic matter. If, in fact, the 
acceleration of the Universe is due to some type of dark 
energy, the big issue that must be answered is: what is 
dark energy? Vacuum energy, some type of scalar matter 
or even some form of more exotic matter? Observational 
data are not enough to allow us decide between the differ¬ 
ent types of dark energy since most of proposed models 
are able to adjust the data seamlessly. We therefore need 
to go deeper in theory in order to achieve a better under¬ 
standing of the mechanisms behind the cosmic accelera¬ 
tion. In this direction, thermal physics is of particular 
importance in dark energy studies. The laws of thermo¬ 
dynamics are based on experimental evidence and apply 
to all types of macroscopic system. Unlike classical me¬ 
chanics or electromagnetism, thermodynamics does not 
predict specific numerical values of observables. Thermo¬ 
dynamics sets limits on physical processes. The power of 
thermodynamics is its generality. Therefore, exploring 
the thermodynamic behavior of the cosmic fluids that 
pervade the Universe may be a line of attack for unveiling 
the nature of the contends of the Universe, particularly 
the hypothetical DE fluid. For example, the positive¬ 
ness of the entropy may be one of the main weapons for 
imposing bounds on the EoS parameter of dark energy 
10]. Applying the laws of thermodynamics to the dark 
energy fluid can help us to constrain or even to rule out 
dark energy models. Below we perform this task. 


I. THERMODYNAMICS OF THE COSMIC 
FLUIDS 

Consider an expanding, homogeneous and isotropic 
Universe filled by matter (baryonic and dark), described 
by a pressure-less perfect fluid (w = 0), radiation, de¬ 
scribed by a perfect fluid with a EoS parameter w = 1/3 
and dark energy, described by a perfect fluid with a EoS 
parameter w = p/pc 2 . Homogeneity and isotropy imply 
that all physical distances scale with the same factor a(t), 
called the scale factor of the Universe. Thus, the physical 
volume of the Universe at a given time is V = a 3 mE 
In such a model the internal energy of the i-th fluid com¬ 
ponent can be written as 

Ui = piC 2 V. ( 1 ) 


1 Here the index 0 will denote the present time value of an observ¬ 

able and we will adopt the convention ao = 1 


Assuming a reversible adiabatic expansion, the first law 
of thermodynamics 

TtdSi = dUi + Pl dV, (2) 

leads to so-called fluid equation, 

dlnpi + 3(1 + Wi)dlna = 0, (3) 

which expresses the energy-momentum conservation. As¬ 
suming that the density is a function of the temperature 
and volume, i. e., pi = pi(Ti , U), the fact that dSi is an 
exact differential implies that 11] 


dlnTi = — Swidhia, (4) 

or, using m to eliminate Wi, 

d In Ti = d In pi + 3d In a. (5) 

Integrating the temperature law ([5]) we get 


Ti_ 



or, in a more suggestive form 


1 PiV 

Wi Ti 


f PifiVg 

Wifi Tifi 


= constant. 


( 6 ) 

( 7 ) 


The above equation generalizes the ideal gas law for a 
time dependent EoS parameter. Finally, the fluid energy 
can be written in terms of the temperature as 


Ui = Uifi^U. 


T 


i, 0 


( 8 ) 


In what follows, we derive the expressions for the heat 
capacity, compressibility and the thermal expansibility. 
These thermodynamical derivatives are easily accessible 
experimentally for any terrestrial fluid. The heat capac¬ 
ity and the compressibility of the fluid are related, re¬ 
spectively, with thermal stability and mechanical stabil¬ 
ity and must be greater than zero. Thus, we can use 
these variables to impose bounds on EoS parameter of 
DE. 


A. The Universe’s heat capacity 

The classical thermodynamical definition of a fluid’s 
heat capacity Ci is 0, 

dQ t = CidTi, (9) 

where dTi is the fluid temperature increase due to an ab¬ 
sorbed heat dQi = TidSi. The heat capacity of a fluid 
will differ depending on whether the fluid is heated at 
constant volume or at constant pressure. From the first 
law of thermodynamics © at constant volume © be¬ 
comes 


dUi = CivdTi, 


( 10 ) 
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where 


B. Compressibility and expansibility 


C iV = 


( dUi 

V 0T t 


(11) If we regard the volume as function of temperatures and 
pressures we have El 


is the fluid’s heat capacity at constant volume. The heat 
capacity at constant pressure can be built up from the 
enthalpy, 



( 20 ) 


hi = Ui+PiV, ( 12 ) 

in terms of which the first law of thermodynamics is writ¬ 
ten as 


dQi = dhi — Vdpi. (13) 

Thus, at constant pressure, m becomes 

dhi = C Pi dTi, (14) 


where 


We define the thermal expansivity, which measures the 
volume thermal expansion at constant pressure, by 


OLi — 


1 

V 



( 21 ) 


and the isothermal compressibility, which measures the 
relative change of volume with increasing pressure at 
fixed temperature, by 


1 (dV\ 

hi dpt) Tt' 


( 22 ) 



is the fluid’s heat capacity at constant pressure. 
From equation ©, it is easy to show that 


Analogously to the isothermal compressibility, we can de¬ 
fine the adiabatic compressibility , if, instead of tem¬ 
perature, the entropy is kept fixed. It can be shown that 
the isothermal compressibility and the isothermal expan¬ 
sibility are related by 


Civ = 777 ^ = constant, 

Ti,o 


(16) 


CX-i 

K Ti 


(dPi\ 

\drjv ’ 


(23) 


for any component of the Universe. Since piV = WiUi , 
the enthalpy m becomes 

hi = (l + Wi)Ui (17) 


and that the ratio between the adiabatic and the isother¬ 
mal compressibilities are equal to the ratio between the 
heat capacities at constant volume and at constant pres¬ 
sure, i.e., 


and, from equations © and (©, we get 

1 din lied 


1 1 , ... Idln^l V, 
C Pi -[l + w i --^—)G iV . 


(18) 


KSi _ CiV 
KTi Cpi 


(24) 


Noting that ptV = WiCiyTi and using © we obtain 


Since t/,;.o = Pi,oc 2 Vo, the specific heat (the heat capacity 
per mass unit) at constant volume is 


OLi — 


C iV / 1 din [wj \ \ 

PiV V 1 3 din a / 


(25) 


C iV = 

Pi, qVq Ti, 0 


(19) 


For relativistic matter T r ^ = 2.725 K so that c r y and 
c Pr are of order of ~ 10 13 cal • g -1 • K _1 . Since the tem¬ 
perature of the others components must be smaller than 
the temperature of relativistic matter, the specific heat 
of the relativistic matter is a bottom limit for the Uni¬ 
verse’s specific heat. As expected, this result reveals that 
the Universe is a huge thermal reservoir. Unfortunately, 
despite easy experimental access for terrestrial fluids, we 
cannot isolate a cosmologically significant portion of the 
Universe, provide an enormous amount of heat and mea¬ 
sure the temperature change of our Universe sample to 
get its specific heat experimentally. 


From (l23ll is easy to show that 


= 


aiV 
WiCiv ’ 


(26) 


and from the above equation and fl^H) we have 


Q^jV 

'UJiCp i 


(27) 


2 Remember that we are assuming that the fluids evolve separately, 
that is, they do not exchange heat, as shown by eq. eg 
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II. CONSTRAINTS ON DARK FLUID 


Thermal stability requires that Civ , C Pi > 0 and me¬ 
chanical stability requires that Ks t , Kt» > 0. Addition¬ 
ally, Civ,C Pi , KSi, andft^ are related by C Pi = + 

j'Yc^ T\^ cx.^ 

—— L and UTi = ns, H— r— L so that C Pj > Civ and 

KTi > KSi- From m and (1251) is easy to see that the 
above conditions are satisfied only if the fluid obeys the 
constraint 


Wi - 


1 din \wi\ 
3 d In a 


> 0. 


(28) 


Reference w w — 1 >0 Va € [0, oo) 


[13] 

W() 

(1 — 6 In a) y 

No 

M 

A wa+/ T 

Wf + 1/ T f 1/T 

a t ' +a- L / T 

No 

[15] 

l , l 

Cl -\~CLj. 

WfWi -, ; 

J w^a 1 -\-w f 

No 

[16] 

Wo + Wo (a - a 2 ) 

No 

[17] 

wo + W'o ^+^2 

No 

[18] 

Wo + w'o 

No 


TABLE I. Thermodynamic viability of some DE parametric 
models found in literature.Here, Wg = (dw/da) a =i. 


It is obvious that if Wi is constant, thermal and mechan¬ 
ical stability implies that Wi > 0 which rule out all neg¬ 
ative pressure fluids with a constant EoS parameter in¬ 
cluding vacuum energy. Thus, if the accelerated expan¬ 
sion of the Universe is caused by a dark energy fluid, 
its EoS parameter must be time-dependent. In Table U 
we list some time-dependent guesses of the DE EoS pa¬ 
rameter. Although in excellent agreement with the data, 
these phenomenological models do not satisfy the ther¬ 
modynamical bound (l28l) . 

Now, let w and pde denotes, respectively, the EoS pa¬ 
rameter and the density of the dark energy fluid. Since 
w must evolves with the time, we can use © to rewrite 
the inequality (E51l as 


3 + din pde 
din a 


dln|iz;| 
— d In a 


(29) 


Integrating both sides of the inequality above we get 


| w | < pde,o _ Tde.o 
wo ~ Pde<J 3 Ide ’ 


(30) 


where we have used ©. According to © the tempera¬ 
ture of the fluid will increase (decrease) with the expan¬ 
sion of the Universe if w < 0 (w > 0). The constraint 
reveals that an eternal accelerated expansion can¬ 
not be sustained. If u> < 0 Va > 1, Tde(o —> oo) —> oo 
and |w(a -A oo)| -A 0 which means that the accelerated 
expansion will stop in the distant future. On the other 
hand, if the sign of w change in the course of the expan¬ 
sion, the Universe can enter and leave in a accelerated 
expansion phase depending on how many times the sign 
of w changes but cannot keep an acceleration expansion 
phase forever. This transient behavior imposed by ther¬ 
modynamics is particularly important for formulation of 
the String/M theory since an eternal accelerated expan¬ 
sion implies that a conventional S-matrix cannot be built 

Emu. 

Now, at present time the inequality (l28l) becomes 

w' 0 > 3wg. (31) 


In order to check the compatibility of the observational 
data with the above thermodynamic constraint we follow 
the approach developed in [22] which is one of the less 
model dependent methods to probe the DE EoS time- 
dependence. This approach consists in assuming that 
the DE density admits a Taylor expansion in the range 
(a — e_, a + e+), that is, 

/ \ f~s , dpnE 

PDE(a) = PDE{a) 


1 d 2 pee 

2 da 2 


da 

(a - a) 2 + 


(a — a) + 


(32) 


and then use the conservation equation © as a recur¬ 
rence formula to write the derivatives of pde in terms of 
the derivatives of w, i. e., 


dpr>E _ 3 

da a 


w)p D i 


d 2 p 


DE 


da 2 


3 . 9 .o 3 dw] 

— (l + w) + —(l+w) 2 -- — 
l a z a z a da 


Pde, 


This approach allows us to constraint w and its deriva¬ 
tives at different redshifts simply changing the series ex¬ 
pansion center a. Since for sufficiently small values of e± 
the second order approximation must work reasonably 
well, we will restrict our analysis up to the second order 
expansion of pue around a = ao = 10 . Thus, choosing 
the expansion center at ao, the second order approxima¬ 
tion of the DE density becomes, 

Pde(o) = PDE.ojl + 3(1 + ui 0 )(l - a) + ^[3(1 + w 0 ) 
-3w(, + 9(l + w;o) 2 ](l-a) 2 }. (33) 


3 For exemple, taking the expansion center a, = ao = 1 and €+ = 
1/3 is possible to cover the redshift range 0 < z < 2. 
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FIG. 1. The wo — w'o parametric space. The thermodynam¬ 
ically forbidden region corresponds to points in the space of 
phase for which the inequality eu is not satisfied. The con¬ 
tours are drawn for Ay 2 = 2.30, 6.17 and 11.8. 


Figure ID shows the observational constraints in 1, 2 and 
3 o' on wo and uj'o f° r a specially flat, homogeneous 
and isotropic Universe filled by relativistic matte, non- 
relativistic matter and a dark energy fluid described by 
obtained from 580 supernovae data of Union 2.1 
compilation, the six estimates of the BAO points given 
in Table 3 of Ref. [26j and the 28 measurements of the 
Hubble function H(z ) compiled by Liao et al. [24] and 
Farooq & Ratra [25| (see also reference therein). The 
present value of the Hubble parameter H$ and the mat¬ 
ter density parameter f2 m ,o were marginalized. For this 
data combination the best fit values are wq = —0.96^21 
and Wq = -0.33^-^ with the upper and lower values 
denoting the one parameter la errors. As we can see, a 
large portion of the u>o — w' 0 confidence regions lies in the 
in the unphysical region. This lack of sensitivity of the 
data to the physical constraint w' 0 > 3wq can be inter¬ 
preted as an evidence against the dark energy hypothesis 
since if dark energy is the true cause of the accelerated 
expansion the data would not be in conflict with its phys¬ 
ical properties. However, if dark energy does not exist the 
data would not be forced to feel its physical properties 
and would force the dark energy parameters to converge 
for the values that better approximates of the true piece 
behind the cosmic acceleration regardless of the physical 
bounds that a hypothetical dark energy fluid must obey. 


III. FINAL REMARKS 

What is causing the accelerated expansion of the Uni¬ 
verse? Is it dark energy or is it that Einstein’s general rel¬ 
ativity that does not work at large scales? Physicists have 
worked on both fronts to answer that question. This is 
the old trial and error method. Particularly, on the dark 
energy side a large number of models have been proposed. 
The diversity of scenarios existing in the literature is due 
to the diversity of aspects to be studied beyond the EoS 
parameter time-dependence, such as fast variation, DE 
implications in the early Universe, and so on. Obviously, 
by adopting this line of attack, nobody expects to find the 
true causes of the accelerated expansion of the Universe. 
Instead it is expected to find clues for the right answer 
between the models’ debris. In this article we believe we 
have taken a big step in understanding the current phase 
of accelerated expansion of the Universe. We have stud¬ 
ied the thermodynamical aspects of an expanding, homo¬ 
geneous and isotropic Universe filled by matter (baryonic 
plus dark), relativistic matter (radiation plus neutrinos) 
and a hypothetical dark energy. We have estimated the 
Universe’s specific heat and examined the constraints im¬ 
posed by classical thermodynamics on the dark energy. 
We have shown that the thermal and mechanical stabil¬ 
ity conditions forbid the existence of negative pressure 
fluids with a constant EoS parameter which excludes the 
vacuum energy as a candidate to explain the cosmic ac¬ 
celeration. We also show that the observational data are 
in conflict with the thermodynamic constraints that a 
general dark energy fluid with a time-dependent EoS pa¬ 
rameter must satisfy. This result suggest that adding 
dark energy to the content of the Universe may not be 
the answer to the cosmic acceleration problem. 

We must noting that, although our analysis excludes 
the vacuum energy, this does not represent the end of 
the cosmological constant. A bare geometrical A-term 
remains in the game if interpreted as a constant of the 
nature whose value must be determined by observations. 
However, what happens to the vacuum energy? Is it null? 
If so, why? We know that the vacuum energy has a sig¬ 
nificant role in the quantum world but, should it play 
any significant role in the Universe at large scales? Can 
we add the quantum vacuum energy so naively to clas¬ 
sical general relativity field equations? We know that 
the fluid description works for relativistic and non rela¬ 
tivistic matter but, can we describe the vacuum energy 
simply as a fluid? Perhaps only a quantum theory of 
gravity can provide an answer to these issues. Beyond 
the issues raised above a question still remains: is the 
geometrical cosmological constant the explanation to the 
accelerated expansion? The A-term certainly is the sim¬ 
plest solution but nobody can guarantee that it is the 
true answer. Thus, finding out deviations of the cosmo¬ 
logical term will remain as one of the hottest theoretical 
investigation lines concerning cosmic acceleration. If the 
dark energy is out of the game, approaches such as the 
kinematic method developed in 27| can be a useful tool 
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to search for such deviations. 

Obviously, we can sacrifice the thermodynamical sta¬ 
bility conditions to keep the dark energy hypothesis alive. 
For example, if we relax the mechanical stability condi¬ 
tion, but keep the thermal stability we have, from (l28l) . 
that w > — 1 saving vacuum energy and quintessence. If 
additionally we give up the thermal stability, phantom 
fields (w < —1) are also allowec0. However, we do not 
think that this is a good way to address the problem. 
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